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Simulation technology can play an important role in rocket engine test facility design and development by 
assessing risks, providing analysis of dynamic pressure and thermal loads, identifying failure modes and predicting 
anomalous behavior of critical systems. This is especially true for facilities such as the proposed A-3 facility at 
NASA SSC because of a challenging operating envelope linked to variable throttle conditions at relatively low 
chamber pressures. Design Support of the feasibility of operating conditions and procedures is critical in such cases 
due to the possibility of startup/shutdown transients, moving shock structures, unsteady shock-boundary layer 
interactions and engine and diffuser unstart modes that can result in catastrophic failure. Analyses of such systems is 
difficult due to resolution requirements needed to accurately capture moving shock structures, shock-boundary layer 
interactions, two-phase flow regimes and engine unstart modes. In a companion paper, we will demonstrate with the 
use of CFD, steady analyses advanced capability to evaluate supersonic diffuser and steam ejector performance in 
the sub-scale A-3 facility. In this paper we will address transient issues with the operation of the facility especially 
at startup and shutdown, and assess risks related to afterburning due to the interaction of a fuel rich plume with 
oxygen that is a by-product of the steam ejectors. The primary areas that will be addressed in this paper are: (1) 
analyses of unstart modes due to flow transients especially during startup/ignition, (2) engine safety during the 
shutdown process (3) interaction of steam ejectors with the primary plume i.e. flow transients as well as probability 
of afterburning. In this abstract we discuss unsteady analyses of the engine shutdown process. However, the final 
paper will include analyses of a staged startup, drawdown of the engine test cell pressure, and risk assessment of 
potential afterburning in the facility. 

Unsteady simulations have been carried out to study the engine shutdown process in the facility and understand 
the physics behind the interactions between the steam ejectors, the test cell and the supersonic diffuser. As a first 
approximation, to understand the dominant unsteady mechanisms in the engine test cell and the supersonic diffuser, 
the turning duct in the facility was removed. As the engine loses power a rarefaction wave travels downstream that 
disrupts the shock cell structure in the supersonic diffuser. Flow from the test cell is seen to expand into the 
supersonic diffuser section and re-pressurizes the area around the nozzle along with a upstream traveling 
compression wave that emanates from near the first stage ejectors. Flow from the first stage ejector expands to the 
center of the duct and a new shock train is formed between the first and second stage ejectors. Both stage ejectors 
keep the facility pressurized and prevent any large amplitude pressure fluctuations from affecting the engine nozzle. 
The resultant pressure loads the nozzle experiences in the shutdown process are small. 

Details of the unsteady shutdown simulation are shown in Figure 1. The shutdown process is initiated by turning 
the mass flow from the engine nozzle inflow boundary off impulsively. The pressure at the boundary is computed 
from solving the momentum conservation equations and does not instantaneously drop. Rather it drops from a 
driving pressure of 1337 psia to 0.16 psia in 12.89 milli-seconds. The evolution of the flow is shown in Figure 1 
where a series of snapshots of Mach number distribution in the facility are shown. As the engine loses power, and 
consequently the pressure drops inside and at the exit plane of the nozzle the strong under-expanded jet through the 
nozzle exit, weakens turns parallel to the jet axis (0.00684 secs). This triggers an expansion wave that travels 
downstream through the supersonic diffuser section. As the expansion wave travels downstream, the strong pressure 
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gradient between the test chamber and the section of the supersonic diffuser close to the nozzle exit drives the flow 
from the test chamber into the supersonic diffuser towards the central axis. The flow entrainment from the test cell 
re-pressurizes the supersonic diffuser section just downstream of the nozzle exit and a small amount of backflow 
enters the nozzle. 

As the rarefaction wave travels downstream near the plane of the first stage steam ejector it changes the shock 
cell structure in the duct between the first stage and second stage steam ejectors. As the incoming flow Mach 
number weakens, it changes the deviation angle of the shock cell that gives the impression that the shock cell 
structure is moving upstream as it is weakening. Flow from the first stage ejectors expands towards the centerline 
and forms a new shock train in the duct between the first and second stage ejectors. Furthermore a compression 
wave is seen traveling upstream from the first stage ejector plane towards the engine nozzle. This compression 
wave first interacts with the expansion wave formed near the corner of the supersonic diffuser where the area 
contraction straightens out. Moving upstream it then interacts with the oblique shock near the diffuser wall where 
the diffuser area contracts. Finally it interacts with the flow expansion between the test cell and the supersonic 
diffuser near the exit plane of the nozzle and remnants of flow from the nozzle forming multiple reflections and a 
complex flow pattern that includes a triple point. 

Meanwhile, as the flow in the duct between the first and second stage ejectors settles into a new shock train 
pattern, the pressurization across the plane of the second stage ejector changes near the centerline in the core flow. 
As a consequence, the oblique shock near the centerline close to the plane of the second stage ejector steepens till it 
turns into a normal shock. Thus the shock train in the duct between the first stage and second stage ejectors 
transitions into a normal shock by the second stage ejectors. More importantly, the flow in the facility past the first 
stage ejector plane remains supersonic throughout the shutdown process. 

Figure 2 shows the line plot of the pressure distribution along the centerline matched against the location and 
pressure distribution in the facility. As a result, the pressure in the nozzle, expansion of the engine plume, shock 
train structure, and pressure drop in exit diffuser can be clearly seen and compared to the physical locations in the 
facility. The line plots (Figure 2 and 3) are plotted in a log scale to accommodate the large variation in pressure 
between the nozzle at full power to the low pressures seen near the nozzle after shutdown. The line plots in Figure 3 
clearly show the pressure dropping in the engine and the downstream traveling rarefaction wave indicated by the 
weakening (smaller amplitudes) of the oblique shocks in the diffuser section just aft of the nozzle. Furthermore, we 
also see the weakening of the shock cell structure in the duct between the first and second stage ejectors and the 
reformation of a new shock cell structure in the duct as the flow from the first stage ejector expands out. The 
formation and eventual migration of the upstream traveling compression wave from the first stage ejectors can also 
be seen as the centerline pressure in the supersonic diffuser rises, as is the transition of the oblique shock near the 
second stage ejectors to a normal shock. Both the first and second stage ejectors are seen to pressurize the duct and 
maintain fairly constant pressure levels. As a result, the pressure distribution in the exit diffuser does not change 
much during the shutdown process. This is indicative of the fact that both the ejector stages are performing 
satisfactorily and prevent any blow-back or large scale pressure fluctuations from entering the engine. Lastly, the 
pressure loading in and around the nozzle after shutdown is fairly small and average pressures range between 100- 
200 Pa (0.015-0.029 psia) as seen in the centerline pressure plots. 
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First stage ejectors 


Second stage ejectors 



Under expanded jet weakens Time = 0.00261 



Rarefaction wave travels downstream Time = 0.00684 



Time = 0.00960 Shock cell structure breaks down 



Flow expands from first stage ejectors Time = 0.01289 



Upstream traveling compression wave Time s 0.01901 Formation of normal shock 



lime = 0.02597 New shock cell structure forms 



Interaction of flow expansion with compression wave Time = 0.03223 



Figure 1. Sequence of Mach Number distributions showing evolution of flow in the facility during engine 

shutdown. 
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Figure 2. Centerline pressure distribution matched up against pressure distribution and physical location in 

the facility. 
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Figure 3. Sequence of centerline pressure plots showing pressure variation along the axis during engine 

shutdown. 
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This paper is the second of a three part series which investigates the computational 
modeling requirements for a new altitude rocket test facility. The focus of this paper is to 
demonstrate the role that simulations can play in supporting the design and development of 
the test facility by assessing risks, providing analysis of dynamic pressure and thermal loads, 
identifying failure modes and predicting the anomalous behavior of critical sub-systems. 
Two separate studies were performed of the facility in different modes of operation. In the 
first study, transient simulations were carried out which represented a nominal engine 
shutdown process. In the second study, the risk of afterburning due to mixing between the 
fuel rich engine plume and oxygen emissions from the chemical steam ejectors was evaluated 
using a chemically reacting simulation. The possibility of facility un-start due to 
afterburning in the diffuser flow path was identified. 


I. Introduction 

The President’s vision for space exploration called for manned missions to the Moon and Mars leading to the 
design and development of two new transport vehicles namely the ARES V and the ARES L The primary 
propulsion element of ARES I and the Earth Departure Stage (EDS) of ARES V will consist of the J-2X engine that 
has evolved from the J-2 engine that was used in the Saturn program. Unlike the SSME engines that are ignited at 
lift-off, flight plans call for the J-2X upper stage engine to be ignited at 133.7 seconds after lift-off at an altitude 
above 177,000 feet. As a consequence, certification and testing of the J-2X engines have to be performed under 
environmental conditions that mimic conditions the engine is exposed to during upper stage ignition and high 
altitude flight operations. There are significant challenges associated with extensive testing of engines of the 300k 
lbf thrust class such as the J-2X engine at high altitude conditions. The A-complex facilities at NASA Stennis Space 
Center(SSC), where SSME altitude tests are performed, are presently limited in scope to altitude testing, and thus a 
new facility (A-3) is under construction to meet the new requirements for altitude testing for the J-2X engine. The 
proposed design includes steam ejectors, a supersonic diffuser, an exit diffuser and a turning duct and is very 
different from existing facilities at SSC. Design issues under consideration with the A-3 facility include: (a) the 
performance of the supersonic diffuser duct and the issues related to “un-start modes” especially under flow 
transients during startup and shutdown, (b) issues related to maintaining low pressure in the test chamber with the 
steam ejector trains, (c) plume aspiration flow dynamics in the test facility that includes flow through the turning 
duct and heat transfer issues related to turning a primarily supersonic plume, (d) entrainment of flow from test 
chamber prior to startup and interaction of initial blast wave with the steam ejectors, (e) purge and pressure control 
in the upstream chamber as well as steam ejector performance especially during shutdown including evaluation of 
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the possibility of blow-back, and (f) performance of chemical steam generators and the interaction of their 
predominantly steam exhaust flow with the test-article plume. 

Simulation technology can play an important role in the rocket engine testing environment by providing design 
support, assessing risks, providing analysis of dynamic pressure and thermal loads, identifying failure modes and 
predicting anomalous behavior of critical systems. In a companion paper 1 , we have demonstrated, with the use of 
CFD, the capability of performing steady analyses of the supersonic diffuser flow path for the purposes of 
evaluating the nominal performance of a sub-scale A-3 facility. However, there are significant transient issues 
related to the performance and operation of the facility especially at startup and shutdown. For example, the 
entrainment of the air in the test chamber by the steam ejectors prior to engine startup is purely a turbulent mixing 
process and is driven by the ratio of the total pressure at the ejector exit to that of the entrained air-stream. Analysis 
based on the ejector driving pressure, ejector nozzle configuration, operating temperature and facility design needs 
to be carried out to evaluate the efficacy of the entrainment process and time required to evacuate the facility to the 
operating altitude pressure levels. On the other hand, the shutdown process in the facility, either as a controlled 
event or a result of an aborted test can result in significant blow-back and pressure loading on the engine resulting in 
serious damage. The placement of the steam ejectors and their operating envelope is usually designed for such cases 
in pressurizing the facility and alleviating any large scale pressure fluctuations that might result from an abort 
related shutdown. 

Lastly, the performance of the steam ejectors/chemical steam generators is critical not only during 
startup/shutdown of the rocket engine but also during the duration of the test. The interaction of the supersonic 
streams, i.e., rocket engine plume with flow from the steam ejector is replete with important flow physical 
phenomena such as shock boundary layer interaction, annular supersonic mixing, shock-shock interaction and 
interaction between the variable area ducting with the supersonic streams. Total pressure ratio between the streams, 
relative Mach number, and the amount of area contraction and expansion play a significant role in determining the 
resultant flow in the facility. Issues with these chemical steam generator units relate to the efficiency of the 
combustion process, resultant emissions and the interaction of any resultant un-bumed oxygen with the fuel rich 
plume in a confined facility. 

In this paper, we will carry out two separate studies with the help of numerical simulations to address transient 
issues with the operation of the facility especially at facility shutdown, and assess risks related to afterburning due to 
the interaction of a fuel rich plume with oxygen that is a by-product of the steam ejectors. The first numerical study 
looked at the engine shutdown event with special emphasis on engine loading. This study was carried out on a 
simplified model of the full scale A-3 facility (Figure la) by removing the turning duct. The second study assessed 
the risk of afterburning in the facility with reacting CFD simulations and for this simulation the sub-scale A-3 test 
facility (Figure lb) was utilized since a facility un-start event was observed during testing in the sub-scale facility. 
In the next section, we provide a brief review of the numerical formulation and the computational framework. This 
is followed by a discussion of the two above-mentioned studies. The paper concludes with a summary of the work 
that was performed and the lessons learned thereof. 




2 nd Stage 
Ejector 


1 st Stage 
Ejector 


Scale 

Rocket 


| Test Cell 



(b) 


Subsonic 

Diffuser 



Figure 1. (a) Schematic of the A-3 Test Stand and (b) the Sub-scale A-3 Test Facility. 
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II. Numerical Framework 


A. Equations and Basic Numerics 

The CRUNCH CFD® code 2 " 3 4 is a hybrid-element (i.e. tetrahedral, prismatic, pyramid, and hexahedral cells) 
unstructured flow solver for viscous, real gas systems. It allows for generalized thermo-chemistry specification, 
permits dynamic grid motion, and has a coupled two-equation turbulence model (see Table I for a list of features). 
CRUNCH is formulated for an edge-based data framework (following Barth 5,6 ) where the solution is saved at the 
cell-vertex and a dual control-volume is defined by cutting across all edges coming to a node. Such an edge-based 
formulation is attractive when dealing with multi- elements since the dual surface associated with an edge can 
include contributions from different element types resulting in a "grid transparent" framework for inviscid flows. 
For efficient computation of large 3D problems, a parallel framework for distributed memory systems has been 
implemented along with a sparse matrix solvers (GMRES and Gauss-Seidel procedures) thereby providing a robust 
computational tool for complex three-dimensional applications. 


Table I. CRUNCH CFD® Navier-Stokes Code Features. 


NUMERICS 

• Finite- Volume Roe/TVD Flux Construction, Vertex Storage 

INTEGRATION 

• Explicit Four-Step Runge-Kutta, Implicit GMRES, Gauss-Seidel 

GRID ELEMENTS 

• Tetrahedral, Hexahedral, Prismatic, Pyramid 

PARALLEL PROCESSING 
CAPABILITIES 

• Domain Decomposition MPI, Independent Grids with Noncontiguous 
Interfacing, Automated Load Balancing 

DYNAMIC GRID 
CAPABILITIES 

• Node Movement Solver (Implicit Elasticity Approach), Automated Embedding, 
Sliding Interfaces 

GRID ADAPTION 

• Variable Element Grid Refinement using Delaunay and cell subdivision 
Procedures, Automated Load Balancing of Adapted Grid 

THERMOCHEMISTR Y 

• Multi-component Real Gas Mixtures, Finite-Rate Kinetics 

TURBULENCE 

RANS/LES 

• ks /EASM Formulations with Compressibility/Vortical Upgrades 

• LES Subgrid Scale Models - Algebraic and One-equation 

• Algebraic (Smagorinsky) and Single Equation (k) SGS Models 

MULTIPHASE FLOW 

• Non-equilibrium Particle/Droplet Solvers (Eulerian and Lagrangian 
Formulations) 

NUMERICS 

• Finite- Volume Roe/TVD Flux Construction, Vertex Storage 

INTEGRATION 

• Explicit Four-Step Runge-Kutta, Implicit GMRES, Gauss-Seidel 


The governing equations are written in finite-volume form for each vertex dual as follows: 

P 

— \QdV+ \F(Q,n)ds = \G(Q,n)ds + \ DdV 

a Q dCi Q ( 1 ) 

Following the standard notation, Q is the vector of dependent variables, F (Q, n) is the inviscid flux vector, G(Q,n) is 
the viscous flux vector, and D is the chemical/turbulent source term. The vectors Q,F,G,D are defined as 

Q= [p,pu,pv,pw,E,pi,...,p s ,...,p NS -i,pk,ps] T (2) 
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(4) 


The first five equations represent global continuity, momentum, and energy equation; the next (NS-1) equations 
represent species continuity; and, the last two equations represent scalar transport of turbulent kinetic energy and 
dissipation rate. 

The inviscid flux procedure involves looping over the edge list and computing the flux at the dual face area 
bisecting the edge. A Riemann problem is solved for using higher order reconstructed values at the dual face (see 
Ref. [2] for details). Presently, a second-order linear reconstruction procedure (following Barth 5 ) is employed to 
obtain a higher order scheme. The higher-order variables need to be limited to yield a TVD scheme. We note that 
the inviscid flux procedure as outlined here is grid transparent; given a dual area associated with each edge the 
details of the different element types contributing to this edge are not relevant. 

The viscous fluxes are computed by estimating the gradients (or stresses) at the cell faces and then performing a 
Green-Gauss integration of the stresses around each dual volume (see Ref. [3] for details). For hexahedral cells 
where the edge vectors are not skewed relative to the cell faces, an edge-based viscous flux procedure can be 
derived. However, for tetrahedral cells where edge skewness is substantial, the edge based procedure fails, 
particularly for the k-c equations, and is reflected in incorrect turbulent viscosity levels. For tetrahedral cells, a cell 
based procedure has been implemented which remedies this problem. 

B. Chemical Kinetics Extensions 

The scalar equation extensions described in the previous section permit the analysis of generalized, multi- 
component gaseous flow-fields including real gas behavior. Source terms added to these equations permit the 
analysis of non- equilibrium phenomena associated with finite rate chemical kinetics. The chemical kinetics are 
modeled in a completely general fashion 2 where an arbitrary set of reactions may be specified by the user. The 
source term vector, D, is written as: 


D = (O, 0, 0, 0, 0, co y cb 2 > 


co 



( 5 ) 
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where co s is the chemical source term for the production of species, s, resulting from any number of chemical 

reactions involving that species. We presently restrict our attention to conventional ’’laminar chemistry” given by 
the law of mass action 


00c 




K f TIa v -K b na v 


( 6 ) 


where a is the molar concentration of species s given by 

= Ps / M s 

The forward and backward reaction rates ( Kf and K^) are a function of temperature and are expressed in modified 
Arrhenius form as: 


K(T)= Qxip(£n C\ +C 2 /T+C 3 £n f) ^ 

where Cl, C2, and C3, are reaction rate coefficients for each reaction. 

The species scalars currently are solved strongly coupled with the fluid equations. For implicit sparse solvers, 
the large block size of the matrices makes this a very memory intensive procedure. A new solution procedure where 
the species convection terms are weakly coupled is currently being investigated. 


C. Turbulence Modeling 

The ‘standard’ high Reynolds number form of the k-s equations forms the basis for turbulence modeling in 
CRUNCH CFD. Transport equations for the turbulent kinetic energy and its dissipation rate are solved along with 
the basic momentum and energy equations. These equations, with supplemental low Re and terms, are as follows, 
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( 9 ) 


where, cr k , o s , Ci and C 2 are the modeling constants, and fi,f 2 ,fju are low Re (near- wall) empirical modeling functions 
which equal unity in the high Reynolds number form. 

Corrections to k-c may be required to account for the effects of compressibility. A compressible wall function 
serves as a baseline method for near- wall analysis but does not perform adequately for highly separated regions. 
Extended variants of k-s that account for non-linear effects may be required to handle the complex vortical physics 
seen in scramjets. Three primary upgrades have been implemented in CRUNCH namely: (1) n + based near wall low 
Reynolds number models to enable integration to the walls; (2) compressibility corrections of Sarkar to account for 
effects of compressible-dissipation and pressure-dilation; and (3) an Explicit Algebraic Stress Model (EASM) to 
include non-linear effects. These extensions are described in reference [4]. The EASM model, in particular, is 
shown to significantly improve predictions in complex vortical flow regions, by allowing the model parameter Cp to 
vary as a function of the local flow conditions, as opposed to being a constant as in the traditional k-e models. 
Developmental work includes accounting for variable Pr t /Sc t effects (see Ref. [7]) and using PDF turbulent 
combustion models (see Ref. [8]). 
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III. Results and Discussion 


A. Simulation of Engine Shutdown in Approximated A-3 Facility 

Unsteady simulations have been carried out to study the engine-shutdown process in the facility and understand 
the physics behind the interactions between the steam ejectors, the test cell and the supersonic diffuser. As a first 
approximation, the turning duct in the facility was removed from the model. The main components of the facility 
such as the test cell, engine/nozzle profile, the supersonic diffuser, the first and second stage ejectors and the exit 
diffuser have been retained and in general correspond to the actual facility with the sole exception of the turning 
duct (see Figure 2). This simplification of the facility was primarily done to facilitate an axisymmetric 
approximation thereby reducing the computational time required and still capture most of the important flow 
dynamics that affect the engine-shutdown process. 



Figure 2. Simplified Model of the A-3 Facility Utilized for Engine-Shutdown Simulation 


The following section of the paper provides a detailed discussion of the flow dynamics observed in the engine- 
shutdown analysis. The engine shutdown process was initiated by impulsively turning off the mass flow from the 
engine nozzle inflow boundary. The pressure at the boundary was computed from solving the momentum 
conservation equations and did not instantaneously drop. Rather it dropped from the nominal combustion chamber 
pressure of 1337 psia to the test cell pressure of 0.16 psia in approximately 12.89 milli-seconds. The evolution of the 
flow is shown in Figure 3 where a series of snapshots of Mach number distributions in the facility are shown. As the 
engine lost its power, and consequently the chamber pressure dropped, the strong under-expanded jet from the 
engine nozzle weakened and turned parallel to the jet axis after 6.84 msec. This triggered an expansion wave that 
traveled downstream through the supersonic diffuser section. As the expansion wave traveled downstream, the 
strong pressure gradient between the test cell and the section of the supersonic diffuser close to the nozzle exit drove 
the flow from the test cell into the supersonic diffuser towards the central axis. The flow entrainment from the test 
cell re-pressurized the supersonic diffuser section just downstream of the nozzle exit and a small amount of 
backflow entered the nozzle after approximately 9.6 msec. 
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First stage ejectors 


Second stage ejectors 



Under expanded jet weakens Time = 0.00261 


Rarefaction wave travels downstream 



Flow expands from first stage ejectors Time ~ 0.01289 



Upstream traveling compression wave Time = 0.01901 Formation of normal shock 



Time = 0.02597 New shock cell structure forms 



Interaction of flow expansion with compression wave Time = 0.03223 



Figure 3. Sequence of Mach Number Distributions Showing Evolution of Flow in the Facility During Engine 

Shutdown 

As the rarefaction wave traveled downstream of the first stage steam ejectors, it changed the shock cell structures 
in the region of the diffuser duct between the first stage and second stage steam ejectors. Effectively, the deviation 
angle of the shock cells steepened which gave the impression that the shock cell structures were moving upstream as 
it was weakening. Flow from the first stage ejectors responded by expanding towards the diffuser centerline and 
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eventually formed a new shock train in the diffuser duct between the first and second stage ejectors. This can be 
seen in Figure 3 at approximately 19msec into engine shutdown. 

The complete expansion of the first stage ejectors to the diffuser duct centerline produced a backward 
propagating compression wave which traveled upstream from the first stage ejector plane towards the engine nozzle. 
This compression wave first interacted with the expansion wave formed near the comer of the supersonic diffuser 
where the area contraction straightens out. Moving upstream it then interacted with the oblique shock near the 
diffuser wall where the diffuser area contracts. Finally, the upstream propagating compression wave had a complex 
interaction with the flow expanding from the test cell near the exit plane of the nozzle and remnants of flow from the 
nozzle at approximately 40msec into engine shutdown. 

Meanwhile, as the flow in the duct between the first and second stage ejectors settled into a new shock train 
pattern, the pressurization across the plane of the second stage ejector changed near the centerline in the core flow. 
As a consequence, the oblique shock near the centerline close to the plane of the second stage ejector steepened till it 
turned into a normal shock. Thus the shock train in the duct between the first stage and second stage ejectors 
transitioned into a normal shock by the second stage ejectors. It is important to note here that in the evaluation of 
these transient flow processes during engine-shutdown, it was evident that the flow in the facility past the first stage 
ejector plane remained supersonic throughout the shutdown process. This signified successful isolation of the 
engine from the external surroundings. If the flow in the diffuser duct upstream of the second stage ejector had 
turned subsonic, then the ambient sea-level pressure of the surroundings would have been able to prematurely 
communicate back to the engine and test cell resulting in a destmctive high pressure loading on the engine nozzle. 

Figures 4 and 5 show line plots of pressure distribution along the centerline of the diffuser duct matched against 
the location and pressure distribution in the facility. The line plots are plotted in a log scale to accommodate the 
large variation in pressure between the nozzle at full power to the low pressures seen near the nozzle after shutdown. 
The pressure line plots clearly show the pressure dropping in the engine and the downstream traveling rarefaction 
wave indicated by the weakening (smaller amplitudes) of the oblique shocks in the diffuser section just aft of the 
nozzle. Furthermore, we also see the weakening of the shock cell structure in the diffuser duct between the first and 
second stage ejectors, and the reformation of a new shock cell structure in the duct as the flow from the first stage 
ejector expanded radially outward to the diffuser duct centerline. The formation and eventual migration of the 
upstream traveling compression wave from the first stage ejectors can also be seen in these figures. Both the first 
and second stage ejectors were seen to pressurize the diffuser duct and maintain fairly constant pressure levels. As a 
result, the pressure distribution in the exit diffuser did not change much during the shutdown process. This was 
indicative of the fact that both the ejector stages were performing satisfactorily and prevented any blow-back or 
large scale pressure fluctuations from entering the engine. Specifically, the pressure loading in and around the 
nozzle after engine shutdown was fairly small, where average pressures ranged between 100-200 Pa (0.015-0.029 
psia) as seen in the centerline pressure plots. 
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Figure 4. Diffuser Centerline Pressure Traces and Corresponding Facility Pressure Distributions at 0.15msec 

After Engine-Shutdown 
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Figure 5. Sequence of Diffuser Centerline Pressure Traces Showing Pressure Variations Along the Diffuser 

Axis During Engine-Shutdown 


B. Afterburning Simulation in Sub-scale A-3 Diffuser 

Reacting flow simulations were carried out for the A-3 subscale facility to study the possibility of afterburning in 
the facility diffuser duct given a fuel-rich engine plume and excess un-bumed oxygen being emitted as a by-product 
from the chemical steam ejectors. Finite-rate chemical kinetics for the H 2 /O 2 system is modeled with a JANNAF 
mechanism comprising of 7 species in 8 elementary reactions. Based on equilibrium calculations using the NASA 
Chemical Equilibrium Analysis (CEA) code, there are 6 primary species (H 2 0, 0 2 , H 2 , H, O, and OH) present in the 
subscale J-2X engine combustion chamber with some additional trace species corresponding to the nominal 100% 
rocket power level operating conditions. In the CFD analysis, the fluid was modeled as a composition of these 
species with the remaining trace species being approximated as an inert species (N 2 ). The CEA code predicted a 
species composition in the sub-scale engine combustion chamber of 90.87% water vapor, 4.773% hydrogen, 3.6% 
OH radicals, and trace amounts of oxygen, atomic oxygen and nitrogen. The above-mentioned species composition, 
combined with a total pressure of 1337psia and stagnation temperature of 5714 deg. F, were the imposed inflow 
boundary conditions in the engine for the CFD reacting flow simulation. The species composition at the exit of the 
converging-diverging steam ejector nozzles for the first and second stages were assumed to be composed of 
primarily water with a non-negligible percentage of excess oxygen. It is important to note here that these simulations 
do not represent the full-scale facility as the chemical steam generators for the full-scale operate using a hydro- 
carbon fuel. Rather, these species concentrations are representative of one of the many tests performed during the 
sub-scale risk mitigation program at NASA Stennis. 
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The mechanism for afterburning in the sub-scale diffuser duct can be explained with the help of Figures 6 and 7. 
Figure 6 shows that in a non-reacting multi-species simulation, the distribution of hydrogen species in the facility 
diffuser duct was primarily concentrated around the central axis of the diffuser until it approached the turning duct. 
Likewise, the oxygen species distribution in Figure 7 shows that un-bumed oxygen exhausted from the first stage 
chemical steam ejector train had its highest concentration close to the walls of the facility. There was as such very 
little mixing between the two supersonic streams (plume from engine and first stage ejector exhaust) in the vertical 
section of the facility. As the flow approached the turning duct, the supersonic core flow in the diffuser (primarily 
hydrogen gas) resisted turning and impinged on the outer wall of the turning duct. The flow separated from the 
inside turn of the diffuser thereby enhancing the mixing between the hydrogen and oxygen. Localized high 
temperature regions along the outer wall of the turning duct, as a consequence of shock formation at the mitered 
joints, acted as ignition sources for the combustible mixture of un-bumed oxygen and hydrogen. A second 
combustion zone was also formed by the Mach disk just downstream of the second stage ejectors. Figure 8 is an 
illustration of the flame fronts that are formed in the turning duct and near the second stage ejectors. These flame 
fronts are illustrated in the figure by iso-surfaces of OH concentration. 



Figure 6. Hydrogen Species Mass Fraction 
Distribution in the Sub-scale Diffuser 


Figure 7. Oxygen Species Mass Fraction Distribution 
in the Sub-scale Diffuser 



Figure 8. Iso-surface of OH Species Concentration Showing the Formation of Combustion (Heat-Release) 

Zones in the Sub-scale Diffuser 
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The OH species distribution along the central plane of symmetry in the facility is shown in Figure 9. There are a 
few important points the authors would like to note first regarding the OH concentrations observed upstream of the 
turning duct. First, given the high temperature field in the core of the engine plume, certain amounts of dissociation 
of hydrogen and water vapor occurred leading to the formation of OH. Secondly, there was a slight increase in OH 
concentration near the shock structures in the supersonic diffuser upstream of the first stage ejectors due to the shock 
induced increase in temperature. Lastly, the high temperature boundary layers also showed formation of OH due to 
dissociation of water. However, none of these OH formation zones represented significant heat release or 
afterburning. 

Although the increase in OH radicals upstream of the turning duct had little effect on the flow-field in this 
region, these radicals did serve to enhance the combustion in the well-mixed oxygen-hydrogen zone within the 
turning duct. The Mach number distribution along the central plane of symmetry is shown in Figure 10. The Mach 
number distribution was plotted on two different scales to highlight the fact that the flow in the turning duct and 
downstream of the second stage ejectors had turned predominantly subsonic leading to un-start of the facility. 



Figure 9. OH Mass Fraction Distribution in the Sub- Figure 10. Mach Number Distribution in the Sub- 
scale Diffuser as Predicted by the Afterburning scale Diffuser as Predicted by the Afterburning 

Simulation Simulation 

A comparison of the temperature distribution (Figure 11) between the afterburning simulation and the 
corresponding frozen simulation (without chemical reactions) shows the impact of heat release due to combustion. 
The temperature in the turning duct has risen considerably (on average by 500K). Furthermore, a comparison of the 
pressure distribution (Figure 12) also shows an increase in pressure in the turning duct (4 psia on average) as well as 
the horizontal duct (6 psia on average). These simulations show that with the present assumed excess concentration 
of un-bumed oxygen emitted from the chemical steam ejectors, there is a potential for un-start in the facility due to 
heat release from combustion in the turning duct. Additional simulations, not reported in the current work, have 
demonstrated that ejector flows with a tenth of the present oxygen concentration produced a much smaller 
combustion zone due to the insufficient oxygen to combust with the excess hydrogen from the engine. Furthermore, 
the low-levels of heat release generated for this condition of low-oxygen concentration steam flow, did not lead to a 
facility un-start condition. 
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Figure 11. Comparison of Temperature Distributions for the Sub-scale Diffuser With and Without 

Afterburning 



Pressure (without afterburning) Pressure (with afterburning) 

Figure 12. Comparison of Pressure Distributions for the Sub-scale Diffuser With and Without Afterburning 

IV. Summary 

This paper is the second paper in a three part series that looked into numerical analyses as a design support tool 
for a steam driven high altitude facility being constructed at NASA SSC. In the first paper, the groundwork was laid 
in terms of modeling approximations for the steam ejector nozzles and the steady-state nominal operating conditions 
of the altitude facility. In this paper, the authors address unsteady issues and assess the risk for facility un-start from 
afterburning. 

Unsteady simulations were carried out to study the engine shutdown process in a simplified facility that was 
devoid of the turning duct. The purpose behind the simulation was to understand the physics behind the interactions 

12 

American Institute of Aeronautics and Astronautics 


RELEASED - Printed documents mav be obsolete; validate prior to use. 


between the steam ejectors, the test cell and the supersonic diffuser. In general, it was observed that as the engine 
lost power, a rarefaction wave traveled downstream that disrupted the shock cell structure in the supersonic diffuser. 
Flow from the test cell was seen to respond by expanding into the supersonic diffuser section and re-pressurizing the 
area around the nozzle. An upstream traveling compression wave emanated from near the first stage ejectors due to 
the first stage ejector expanding to the center of the duct. However, both stage ejectors kept the facility pressurized 
and prevented any large amplitude pressure fluctuations from affecting the engine nozzle. Overall, the resultant 
nozzle pressure loads due to the nominal shutdown process were small and within the acceptable range. Of course, 
the pressure waves predicted during engine-shutdown were unique to this model, especially considering the fact that 
the turning duct was not included. The full diffuser system has been modeled, yet it’s behavior is dependent on the 
GN2 test cell purge rate specified. Further investigation is currently planned to model the full system with the 
expected GN2 rates and ejector shutdown timing as completion of a Phase II SBIR program in this area. 

Reacting flow simulations were also carried out in the current work to study the risk of afterburning in the 
facility and provide any necessary recommendations for mitigating potential problems. The numerical simulations 
revealed that in the case of excess unbumed oxygen being emitted from the steam ejector units, there existed the 
possibility of afterburning inside the diffuser duct. The excess hydrogen from the engine and oxygen from the 
ejectors mixed in the turning duct and underwent shock-induced ignition both along the outer wall of the turning 
duct and at the exit of the turning duct near the second stage ejectors. As a result of the concomitant heat release, 
the pressure and temperature in the turning duct rose resulting in a facility un-start condition. Further 
simulations/investigations are underway to determine a safe threshold for oxygen concentration where the facility 
can operate without combustion facilitated disruption. 
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